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Abstract. - The dynamic behaviour of glassy materials displays strong nonequilibrium effects, 
such as ageing in simple protocols, memory, rejuvenation and Kovacs effects in more elaborated 
experiments. We show that this phenomenology may be easily understood in the context of 
the nonequilibrium critical dynamics of non-disordered systems, the main ingredient being 
the existence of an infinite equilibrium correlation length. As an example, we analytically 
investigate the behaviour of the 2D XY model submitted to temperature protocols similar to 
experiments. This shows that typical glassy effects may be obtained by 'surfing on a critical 
line' without invoking the concept of temperature chaos nor the existence of a hierarchical phase 
space, as opposed to previous theoretical approaches. The relevance of this phenomenological 
approach to glassy dynamics is finally discussed. 

...comme des amateurs de surf qui glissaient sur la vague. Vassili Axionov, Une brulure. 

Experiments on a large class of glassy materials demonstrate intriguing similarities in 
their nonequilibrium dynamic properties. Structural — in particular polymeric — glasses are 
the paradigm, and have been studied in a systematic way over a long period time Interest 
in ageing experiments grew with the realization that very similar effects are observed in spin 
glasses M. More recently, the family of glassy systems has been much enlarged to include 
systems as different as dirty type-II superconductors [0, several complex fluids || , disordered 
dielectrics |], [| and ferromagnets Q , electron glasses || , and granular materials || . 

These systems have in common that, in a part of their phase space, they cannot be equili- 
brated within the experimental time window: they undergo a 'glass' or 'jamming' transition. 
As a result, they exhibit nonequilibrium dynamics, where the whole history of the sample 
becomes relevant The conceptually simplest experiment performed in glassy materials 

consists of a sudden quench of the system into its glassy phase at the initial time t w = 0. The 
resulting ageing behaviour is best observed through the measurement of a two-time quantity, 
C(t,t w ), typically a correlation function or a linear susceptibility. Experiments show that 
time translation invariance is absent, i.e. C(t,t w ) ^ C(t — t w ). Remarkably they also show 
that the generic scaling behaviour C(t,t w ) ~ C(t/t w ) is very often observed, implying that 
the relaxation time t re i of the system is of the order of its age, t Te \ ~ t w . 

This ageing phenomenology can be explained by models of coarsening of domains of equiv- 
alent thermodynamic phases |l^,[Tl]]. The canonical example is a ferromagnet quenched into 
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its ordered phase which exhibits ageing, due to the growth of ferromagnetic domains of dif- 
ferent orientations. A length scale irit), over which equilibrium is established, grows with 
time t at temperature T, leading to the above scaling behaviour [ jnj . However, such simple 
models fail to explain the spectacular 'memory' and 'rejuvenation' effects observed in 
more elaborate experimental protocols unless complemented in some way. One approach is to 
include temperature chaos [13] , which amounts to a great sensitivity of equilibrium states to 



temperature changes, as first postulated in the context of spin glasses 14 . The very existence 
of temperature chaos in this sense is however still debated, with rather negative conclusions. 

Ageing may alternatively be viewed as the search by the system of the deepest 'traps' of 
a complex energy landscape ]l5| . Extending this approach to a hierarchical structure of the 
traps allows for a qualitative understanding of temperature cycling experiments [JL6| , Multi- 
trap models of this kind are studied in Ref. jl7[ . This approach has however the disadvantage 
of not providing a real space description of the problem. Ageing is also captured in a more 
microscopic way by the exact solution of several mean-field disordered models These 
do account for rejuvenation and memory |l8| ], but the explanation relies on the notion of 
dynamic ultrametricity, which is incompatible with experimental data |T§| ]. 

In this paper we show that a phenomenological domain growth approach may be used 
to interpret T-cycling experiments, provided the dynamics is considered at criticality. For a 
quench into an ordered phase, the equilibrium correlation length £ cq is microscopic and there- 
fore decoupled from lr{t)- However, the key point here is that quenching to a critical point, 
one has £ oq = oo, with the result that equilibrium fluctuations occur on all length scales up to 
the dynamic correlation length = Ixit). Temperature cycling along a line of critical points 
thus affects critical fluctuations of all length scales between the lattice spacing and the system 
size. This allows for both rejuvenation and memory effects without making the assumption 
of temperature chaos, or evolution in a hierarchical phase space. This may also be viewed as 
a concrete realization of the multi-length scale scenario discussed by Bouchaud As an 

example, the 2D XY model is analytically considered, although the physical interpretation of 
the results is of more general validity. Finally, we argue that this phenomenological approach 
may be relevant to many glassy materials. 

Ageing at a critical point. - The 2D XY model presents a line of critical points in the 
temperature range T € [0, Tkt] where Tkt is the temperature of the Kosterlitz-Thouless 
transition. The model describes a two-component order parameter ip(x, t) in a bidimensional 
space. The properties of the low temperature phase, < T < Tkt, can be described by the 
spin-wave approximation Jl9| 

H[6] = e@.J d^lV^I 2 , (1) 

with ip = e l8 and where p(T) is the spin-wave stiffness. The dynamics is modelled by a 
Langevin equation 

dO(x,t) _ 5H[6] 

at ~ S9(x,t) +ax > t} - (2) 

The last term is the thermal noise, described by a Gaussian variable of zero mean and variance 
(C(x, t)C(x', t')) — Attt](T)p(T)5(x — x')S(t — t'), where rj(T) is the usual critical exponent, 
linked to the spin-wave stiffness by 2irrj(T)p(T) = T. The dynamic behaviour is analysed 
through the two-point correlation function C(r, t) = (<p(x, t) ■ tp(x + r, i)), and the two-time 
autocorrelation function C(t,t w ) = {ip{x,t) ■ <p(x,t w )). Due to the Gaussian character of 
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the Hamiltonian (|l|), these quantities follow straightforwardly from the computation of the 
angle-angle correlation function Cg(k,t,t w ) — (8(k,t)9(—k,t w )). 

The ageing dynamics of the model was discussed in details in Refs pO|-p^|. In particular, we 
recall that topological defects are not described by the spin wave approximation ([!]) , although 
they can influence the dynamics following a quench pl| . We thus choose initial conditions 
such that no vortices are initially present in the system, with Cg(k, 0, 0) = 2irrj(Ti) / 'k 2 , where 
Ti is the initial temperature This implies C(r,0) ~ {r/a)^ ri ^ Ti ^ at large distances. The 

,2 2 

UV cutoff a is introduced through the factor e a in all integrals over the Fourier space, 
simulating the lattice spacing. Although the correlation length in the initial state is not 
defined, the correlation function decays very rapidly if T}{Tj) is chosen to have a large value. 

The dynamics following a quench to temperature Ti at time t = is solved by Fourier 
transforming (^), using (|l|). For the angle-angle correlation one obtains 

C e (k, t, t w ) = ex P [-(H Tl (f - t w )f] + 2 ^ 2 ~ Vl) exp[-(M Tl (t + t w )) 2 ]. (3) 

We define the dynamic correlation length £r k (t) — (pkt) 1 ^, z = 2, together with the notation 
rjk = ^(Ife) and pk — p(Tk). Two-point and two-time correlation functions easily follow from 
(|h. For times such that It x (t — t w ) » a, one finds 



C(r,t) ~ exp 
C(t, t w ) 



\2V2£ Tl (t) 



iT^t-t^J \4£ Tl (t)£ Tl (t 



(4) 



where G(x) — | + 2 lnx + dt ^— j ; C is the Euler constant. For distances a <C r <C 

£Ti(t), equilibrium behaviour is obtained at Ti, C{r,t) ~ (r/a) _I?1 , while at large distances, 
^Ti (i) -C r, one has C(r, t) ~ (r/a)~ rii {ajtr-i {t)) r,1 ^ r,i , which is the same power law behaviour 
as in the initial state, with a time-dependent multiplicative factor |22|. Similarly, the two-time 
correlation function displays an initial equilibrium behaviour C(t, t w ) ~ (t — t w )~ ,n / 2 , which is 
interrupted after a time t — t w ~ t w , when the relevant scaling variable becomes ^ (t ) /£ti ■ 
The dynamics is thus interpreted as consisting of two types of critical fluctuations pl|.|2~3]]. 
For wavevectors kl^^t) 3> 1, the equilibrium state at T\ has been reached, while fluctuations 
with k£Ti(t) <C 1 arc still very near to their initial noncquilibrium state. The ubiquitous t/t w 
scaling is thus obtained as a consequence of a growing equilibration length in the system, as 
in ordinary coarsening [^o| . The only difference is that here the correlation length is fixed to 
£r{t), while in the former it is microscopic. 

Rejuvenation and memory effects. - In temperature cycling experiments, the system is 
first quenched to a temperature T± where it ages. After a time t\, the temperature is shifted 
to T 2 < Ti. It is observed that the ageing is restarted by the temperature change: this is 
the rejuvenation effect. Then, at time ti > t\, the temperature is changed back to T\. After 
a short transient time the ageing is observed to be the same as if t% = ti, that is, as if the 
temperature cycle had not existed: this is the memory effect (see Fig. [j]). The coexistence of 
these two effects was popularized by the 'dip-experiment' proposed in Ref. p4| , which has since 
been largely used, showing again very similar results in different materials |2|,|^,p|-[7|,[l6|,p5|. 

The rejuvenation effect can be conveniently demonstrated by taking the limit t\ — > oo, in 
which case the system starts from an equilibrated configuration at T± and is quenched to T 2 . 
This is in fact precisely the situation considered in the previous section if the changes Tj — > Tj 
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Fig. 1 - Evolution of the imaginary part of the Fourier transform of the autocorrelation function in 
a temperature cycling. Parameters are: a = 1, tv = 50, r\i = 2.0, t\ — 6000, 771 = 0.7, p\ = 1.0, 
t2 — 12000, r)2 = 0.35, p2 = 0.05. Inset: the time interval at T2 has been dropped out. To better 
demonstrate the memory, C"(u,t) — C"{lu, 00) has been plotted in a lin-log scale, the full line being 
the behaviour with t2 =ti, i.e. without the cycle. 



and Ti — > T2 are made. Thus, Eqs. ([|) illustrate the fact that ageing is restarted when the 
temperature is changed along the line of critical points. Keeping t\ finite, the angle-angle 
correlation becomes, at times t > t w > t\, 

C e (k, t, t w ) J^- exp[-(Wr a (* - O) 2 ] + Mm ~ m) exp[~(ke T . 2 (t + t w - 2^)) 2 ] 

(5) 

+ 27rfa-r ?1 ) exp[ _ (HTa(f + tw _ 2h))2] exp [„2(H Tl (t 1 )) 2 ]. 

When t\ — > 00 , the last term vanishes and the results of the previous section are indeed 
recovered. On the contrary, when lr 2 {t ~ ii) ^ &Ti(ti), the last two terms combine to give 
an ageing behaviour as if the system was directly quenched to T2. Equation (^J) shows that 
fluctuations with ki^ (ti) S> 1 age at Ti as if they were quenched from T\ (the last exponential 
makes the last term vanish), whereas those with fc^Ti(^i) <C 1 age at T2 as in a quench from 
Ti. In this case the last two terms may indeed be combined when the last exponential is ~ 1 
to make the contribution proportional to 7/1 disappear. 

Physically, rejuvenation results from the fact that when the temperature is shifted to 
T 2 <Ti, all fluctuations fall out of equilibrium. Fluctuations such that ki^iti) ^> 1 have to 
adapt their Boltzmann weight to the new temperature, while for fc^Ti(ii) < 1, no equilibrium 
had been found at the previous temperature. Thus even when t\ = 00, ageing is restarted by 
a temperature change. We emphasize that the effect is absent in a standard coarsening below 
T c , where the reequilibration of thermal fluctuations following the shift to T2 is very fast, with 
an associated time scale t r given by £T 2 (tr) ~ £cq(7i) "C ^Ti(^i)- 

The memory effect can be seen as the temperature is changed back to Ti at total time i 2 , 
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with ^t 2 (^2 — ti) <C ^Ti(ii)- The angle-angle correlation function reads now 

Ce ( k , t, t w ) =-^ e -(«n e-*»)) 2 + Mm-m) e - { u Tl (t+t m -2t 2) f l _ e -2 ( ki T2 {t2 -H)f 

k 2 k 

| 27rfa - m) c -2(ki T , (ti)) 2 e -2(M T2 (t 2 -t 1 )) 2 e -(fe£ Tl (t+t„-2t 2 )) 2 ^ 

k 2 

(6) 

To understand this result, three types of fluctuations have to be considered. 

(i) k^ 1 <C ^t 2 (*2 — These fluctuations are equilibrated at T 2 and have to reequilibrate 
at T\. This can be seen by setting e ~ 2 ( M T 2 (* 2 -ti)) ^ o in Eq. from which one finds an 
equation similar to (||) , describing a shift from T2 to T\ . The reequilibration takes place on a 
time scale t m such that t^itm) ~ ^T 2 (h — t\) <C fri(ii), i-e. t m <C t%. This is the short initial 
transient observed in experiments and in Fig. 1. 

(ii) lT 2 (t2 — t\) <C fc _1 <C #ri(ii). Only the first term of Eq. (Q) survives. These fluctuations 
had no time to equilibrate at T2 but were equilibrated at T% before the cycle. Hence, they 
undergo equilibrium dynamics at T± immcdiatly after the cycle. It is in these length scales 
that the memory is conserved. 

(iii) ^Ti(ti) <C fc _1 . Only the last term in Eq. @ survives for these wavevectors. Taking 
e -2(fef Tl (ti)) ^ ^ anc j e -2(fc£ T2 (t 2 -ti)) ^ Eq. is recovered, describing a quench from Tj to 
Ti . These length scales retained their initial nonequilibrium state throughout the temperature 
cycle and will thus equilibrate at Tj through the further growth of It x . 

As a consequence, after the transient time t m , the dynamics appears as if the cycle at T 2 
was absent, accounting for the memory effect observed in experiments. We show in Fig. [I] 
the behaviour of the imaginary part of the Fourier transform of the autocorrelation function 

in a temperature cycle. At equilibrium, this function is directly 
related by the fluctuation-dissipation theorem to the linear susceptibility which is the 

quantity measured experimentally. The figure demonstrates that the rejuvenation and memory 
effects are quantitatively accounted for by the 2D XY model. The two-time autocorrelation 
function is computed from Eqs. (|^) and @, as was done to obtain Eq. (^) from Eq. (||). As 
it is rather complicated, we do not report it here, its interpretation being the same as that 
developed above. 

Interestingly, the same effects could be observed in a positive T-cycle, T% > T\, provided 
^t 2 (*2 — ti) <C £j\ (ti). This feature is in qualitative contradiction with the hierarchical phase 
space picture fTq| , and experiments could thus distinguish between the two approaches. 

Kovacs effect. - In the literature of structural glasses, 'memory effect' has a different 
meaning p|,p6[. Here, the system is first quenched to a given temperature T\. The slow 
nonequilibrium evolution of a one-time quantity V(t) (volume, refraction index) is followed. 
The temperature is then raised at time t\ to a value T 2 such that, at equilibrium, V would 
have the value it has just before the shift, i.e. V(ti) — V cq (T2). Whereas it would be naively 
expected that V(t > t\) = const, a nontrivial non-monotonic behaviour is observed instead. 
To avoid confusion, we will refer to this effect as the Kovacs effect [^6|. Again, we note that 
it was observed in several materials [| pi [ffi-p9j . 

Within the 2D XY model, one follows the behaviour of the energy density e(t) in a protocol 
Ti -> Ti at t = and T x -> T 2 at t x , with T l >T 2 >T 1 . For times such that £ T2 (t - t x ) > a, 
we get in the harmonic limit 

m - _L ( 2 Il 4. T *~ Tl 4. \ (7) 

e[t) igttU 2 4 2 (*- *i) + 4x(*0 e 2 T 2 (t-h)J- {) 
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Fig. 2 - Kovacs effect in the 2D XY model. Full lines are for various intermediate temperature 
7i = 0.95, 0.9, 0.8, 0.7, 0.6, 0.5 (from bottom to top), with T t = 5.0 and T 2 = 1.0. 



The energy first increases due to the third term, since r\\ < r\i- When £r 2 {t — ti) ^ ^Ti(*i)> 
the second and third terms combine to give the simple ageing behaviour of a system directly 
quenched from Tj to T2, and hence a decreasing energy density which merges with the curve 
obtained in this simple experiment. Equation (Q) was used to build the curves in Fig. ||, which 
reproduce quantitatively experimental findings. 

Immediately after the shift at time ti, fluctuations with fc _1 <C £ti(*i) are equilibrated 
at T\ and have to increase their energy in order to equilibrate at the new temperature T2. 
Fluctuations with ^Ti(^i) "C fc _1 , which retained their nonequilibrium state at Ti, decrease 
their energy to reach their equilibrium state at T2. These two types of fluctuations have 
opposite contributions to the energy evolution and act on different time scales. Thus, a 
maximum occurs at a time tk> such that £r 2 (ik — ti) ~ ^Ti(t{), the larger T2 — T\ the larger 
the height of the maximum. This interpretation can be viewed as a spatial transposition of 
the distribution of time scales usually invoked to account for this effect (TJ . 

Discussion. We have shown that coarsening models at criticality reproduce generic 
ageing, rejuvenation, memory and Kovacs effects seen in many glassy systems, thus providing 
a new and simple phenomenological domain growth approach to ageing phenomena. 

This seems to beg the question: are glassy systems critical? A possible answer is yes they 
are, in which case glassy systems should have in common a line of critical points below the 
glass transition, in analogy with the 2D XY model. This is not inconsistent with existing 
theoretical descriptions of the spin glass phase but seems less likely in other systems, 
such as disordered ferromagnets JtJ or ferroelectrics pj. Another is no; but they do at least 
behave as if there is a line of critical points over the experimental time window. The second, 
more pragmatic solution requires that the coherence length irit) is never decoupled from 
£ eq (T), even when the latter is finite. One might argue this if the excitations on length scales 
exceeding £ oq (T) are dominated by activated processes, giving exponentially increasing time 
scales for £r(t) > £ q- These conditions would make it virtually impossible to enter a regime 
where tx{t) 3> t; cq (T). In this case, a slow evolution in the crossover region could result in 
effective critical behaviour with continuous evolution of the exponents, with the temperature, 
which is the main ingredient of our analysis of the rejuvenation effect. This is fully consistent 
with results in 3D spin glasses [|9|,^|. It also implies that similar results could be obtained 
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in spin glass films around the T = critical point |3l]] . 

In both cases, of course, disorder is essential. Indeed, we would not expect to observe 
ageing phenomena easily in a non-disordered critical system as the divergence of time scales is 
actually very weak and can in general be beaten easily, both experimentally and numerically. 
Disorder would therefore be needed to provide the necessary separation of time scales jl5| and 
to put the ageing phenomena discussed here within the observed experimental time window. 
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